Optimization of the leak conductance in the squid giant axon 
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We report on a theoretical study showing that the leak conductance density, Gl, in the squid 
giant axon appears to be optimal for the action potential firing frequency. More precisely, the 
standard assumption that the leak current is composed of chloride ions leads to the result that the 
experimental value for Gl is very close to the optimal value in the Hodgkin-Huxley model which 
minimizes the absolute refractory period of the action potential, thereby maximizing the maximum 
firing frequency under stimulation by sharp, brief input current spikes to one end of the axon. 
The measured value of Gl also appears to be close to optimal for the frequency of repetitive firing 
caused by a constant current input to one end of the axon, especially when temperature variations 
are taken into account. If, by contrast, the leak current is assumed to be composed of separate 
volt age- independent sodium and potassium currents, then these optimizations are not observed. 

PACS numbers: 87.19.Ib,87.19. 11,87.19. lo,02.60.Pn 
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I. INTRODUCTION 

A considerable amount of evidence has emerged in re- 
cent years to show that many of the parameters which 
govern the structure and function of biological ner- 
vous systems are at optimal values for metabolic en- 
ergy consumption, information rates, or some combina- 
tion thereof, presumably because of evolutionary pres- 
sures Many of these studies have focused on the 
squid giant axon because of its well-known and relatively 
simple properties. Hodgkin and Adrian hypothesized as 
early as the 1970s that the channel densities in the squid 
giant axon are at values that maximize the action poten- 
tial velocity [1, Q , although more updated axon models 
have called this into question @ and suggested that axons 
are optimized for the energy associated with the action 
potential instead 0, [1] • There are a large number of in- 
dependently variable parameters that significantly affect 
the functioning of the squid giant axon, only a few of 
which have been systematically investigated for possible 
optimizations. Here, we present results for one of them: 
the leak conductance. 

The voltage-independent leak conductance is one of 
three conductances known to be present in the squid gi- 
ant axon. With a measured value of about 0.3 mS/cm^, it 
is much smaller in magnitude than the maximum voltage- 
gated Na+ (120 mS/cm^) and K+ (36 mS/cm^) conduc- 
tances, yet it nevertheless plays an important role in the 
electrical stability of the axon. Leak conductances are 
known to be present in many other kinds of neurons as 
well, such as molluscan pacemaker cells [1]. 

Because of its small size, there has long been de- 
bate about the exact nature of the leak conductance, 
with some of the debate centering on how much of it 
is due to voltage-gated K"^ and Na"*" channels that re- 
main open at rest [l^, Ull- However, a nonselective. 
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voltage-independent cation channel protein has rece ntly 
been conclusively identified in mammalian neurons [13 |. 
In squid giant axons, the closeness of the leak reversal 
potential to the equilibrium potential of Cl~ has tradi- 
tionally been taken to indicate that chloride ions are a 
significant contributor to the current, although there may 
be others [isj . 

Given the possibility that the leak conductance is an 
easily evolvable parameter capable of influencing the elec- 
trical properties of axons, it is natural to ask whether the 
leak conductance is at an optimal value for some quantity 
related to information processing or energy consumption. 
Here, we present results showing that the leak conduc- 
tance is near-optimal for the absolute and repetitive firing 
frequency of the axon if the leak current is assumed to 
be chloride, but not if it is sodium/potassium. For com- 
pleteness we also investigated the effects of Gl on the 
relative refractory period, which is substantially harder 
to calculate numerically. We did not find strong evidence 
that Gl is optimal in this case for either channel model. 



II. METHODS 
A. The Hodgkin-Huxley Model 

The squid giant axon, about 0.5 mm in diameter, is 
one of the largest axons in nature. It innervates muscles 
in the squid mantle, and its action potentials cause the 
muscles to contract, expelling a brief jet of water and al- 
lowing the squid to move away quickly from danger. The 
axon is postsynaptic to neurons in the dorsal magnocel- 
lular lobe, and it is indirectly connected to the ventral 
magnocellular lobe, which integrates sensory input. 

Fortunately, the squid giant axon is also one of the sim- 
plest known axons, being unmyelinated and having only 
two voltage-gated ion channels with relatively straight- 
forward kinetics. As such, we arguably know more about 
the squid giant axon than any other neural system, and 
it is possible to model it to a high degree of biological 
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accuracy. The Hodgkin-Huxley (HH) model, based on 
the experiments of A. L. Hodgkin and A. F. Huxley in 
1952 M, remains very useful, although subsequent re- 
finements to the sodium and potassium channel kinetics 
have been made [HI, [l^] . In our study, we used the tra- 
ditional version of the HH model described below. 

The model treats the squid giant axon as a cylinder 
of length L and uniform diameter d. For most of our 
simulations, L was set to 0.8 cm. The diameter was gen- 
erally in the range of 300-600 |J.m. (These dimensions 
are typical of the biological ones.) The cylinder has an 
axial resistivity Ra representing the axoplasm. This was 
generally set to 35.4 ■ cm after Hodgkin and Huxley's 
measurement. Transmembrane currents flow through 
two voltage-dependent Na'^ and conductances, Gn.-, 
and Gk, and a voltage-independent leak conductance, 
Gi. The Hodgkin-Huxley experimental measurements of 
Gi ranged from 0.13 to 0.5 mS/cm^, with an average of 
about 0.26 mS/cm^[13]. It was primarily this value that 
we varied in our simulations. The voltage-dependent con- 
ductances are functions of time and, indirectly, the mem- 
brane potential; 



G«., = G«. [mit)fh{t), 



Gk = Gk [n{t)Y 



(1) 



(2) 



where m and h are state variables representing the frac- 
tion of Na^ channel subunits in the open and non- 
inactivated states, respectively, and n is the state variable 
representing the fraction of open channel subunits 
(squid K"*" channels do not inactivate) . Each of the four 
subunits in a channel has to be open or non-inactivated in 
order for the channel to pass ions. Thus, m^h and are, 
respectively, the fraction of Na"*" and K"*" channels which 
are open. These state variables evolve according to equa- 
tion 2] below. The maximal conductances, Gn^ and Gk, 
obtain when all the channels are open. We used Hodgkin 
and Huxley's experimental values of Gn., = 120 mS/cm^ 
and Gk = 36 mS/cm^. 

The cell membrane has a constant intrinsic capacitance 
of approximately Go = 0.88|^F/cm^ [l3]. The voltage- 
gated sodium channels also contribute a phenomenologi- 
cal capacitance, Gq = 0.13 |J.F/cm^, the so-called "gat- 
ing" capacitance |5|. (In principle, the voltage-gated 
potassium channels contribute a gating capacitance too, 
but this is so much smaller than the sodium gating ca- 
pacitance that it can be neglected.) 

Electrical excitations of the axon are described by four 
coupled differential equations. The first of these is a mod- 
ified version of the cable equation: 

+G^n'(V - E^) + G,(V - E,) (3) 

where V is the cross-membrane potential, with the extra- 
cellular side taken as ground. The Na"*" and reversal 



potentials, i?Na and i?K, are determined by the ionic con- 
centration gradients across the membrane: the values we 
used, -Enb = 50 mV and — -77 mV, are typical of 
the squid giant axon. The leak reversal potential, E^^, 
is determined by the ion(s) which pass through the leak 
channels and is experimentally around -55 mV [T3 |. We 
discuss our leak channel assumptions in greater detail 
below. 

The other three differential equations govern the gating 
variables: 



ds 
di 



= as{V) X (I ^ s) - l3siV) X s , 



(4) 



where s — n, m, or h. The rate coefficients, as 
and /3s, are the fraction of s subunits per unit time 
switching from closed/inactivated to open and open to 
closed/inactivated, respectively. These rates were empir- 
ically measured by Hodgkin and Huxley as: 

I3rn iV) = ^x 0.4 X e-(^+65)/i8 [ms-i] (6) 
ah iV) = (f>x 0.07 X e-(^+65)/20 [ms-i] (7) 

^''(^)-^X (,-(V+35)/10 + l) (8) 

aniV) = c^xOMxj-^^}^±§L_ [nis-] (9) 
^„(y) = X 0.125 X e-(^+65)/80 [^g-ij (^0) 
with V in mV. The temperature coefficient (j) is 

^ ^ 3(T-6.3)/10 (11) 

with T in °C. It should be noted that the openings and 
closings of individual ion channels are stochastic in na- 
ture; Eq. ^ describes the average behavior of a large 
ensemble of s-subunits. 

Eqs. ^ and (g]) are highly nonlinear (Eq. (4) be- 
cause of the forms of the rate coefficients (5)-(10)); and, 
without significant approximations, they are analytically 
intractable. However, it is known that they have a unique 
solution describing a single voltage spike, or "action po- 
tential," propagating at a uniform velocity Vap along the 
axon. The action potential velocity is a function of the 
different biophysical parameters in Eqs. (CTIU]). though 
it does not have an exact analytical form and must gen- 
erally be either approximated or determined numerically. 

Our simulated axon contained 1000 isopotential seg- 
ments, each of length 100 |xm. Equations [3] and S] (for 
n, m, and h) were solved simultaneously in each segment 
using an implicit backward Euler method. Our time step 
was 1 \JiS. We verified that the time and spatial resolu- 
tions were sufficiently fine so as to not significantly infiu- 
ence our results. We assumed T = 18.5 °C unless stated 
otherwise. 
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B. The Leak Channel 

We tested two different assumptions about tiie nature 
of the leak channel. We assumed first that it is a voltage- 
independent Cl~ conductance, in which case the leak re- 
versal potential, i^i, is just the equilibrium potential of 
chloride: 



E,, = Er 



-55 mV 



(12) 



We identify this case as "Cl~ leak" in the figures. In 
our simulations, we typically varied the value of Gi while 
keeping all the other parameter values the same. This has 
the effect of changing the resting potential of the axon 
(see Figure [l} and, indirectly, the maximum frequency at 
which action potentials can fire. We numerically calcu- 
lated the new resting potential and set the axon to this 
value at the beginning of our simulations before action 
potentials were evoked. 

In the second case, we assumed that the leak channel 
consists of two voltage-independent Na"*" and K"'" con- 
ductances, in which case Ei^ is determined by a weighted 
average of E^^ and E^- 



E,= 



(13) 



In this case, which is identified as "Na~'~/K''' leak" in the 
figures, the total leak conductance is just the sum of the 
sodium and potassium leak conductances: 



G, — G, 



G, 



(14) 



Since only Gi and Ei, appear in the equation of motion 
for V, Eq. ([3]), the refractory periods are not affected by 
whether it is Cl~ or Na^/K^ that goes through the leak 
channel if only G^ is varied. That is, any optimization 
results involving the leak conductance alone that obtain 
for Cl~ leak channels should obtain for Na~''/K"^ leak 
channels as well - or, in fact, for a leak current composed 
of any combination of permeant ions, provided that their 
overall reversal potential is approximately Eqi- In order 
to more completely distinguish between the two cases, 
we added the further requirement in the Na^/K^ leak 
channel case that whenever the value of G^ was altered, 
the value of Ej^ was altered as well so as to keep the overall 
resting potential, Er, at -65 mV. This is equivalent to 
altering the ratio of the Na"*" leak conductance, G^Na, to 
the K"'" leak conductance, GiK, while keeping the reversal 
potentials E^^^ and -Ek the same. 

Thus, the mathematical distinction between the two 
models is that for CP leak channels, only G^ is varied 
in Eq. (j3|), which in turn causes the resting potential to 
vary as shown in Fig. ([Ij. For Na^/K"*" channels, by 
constrast, both Gi, and E^, are varied in tandem so as to 
keep the resting potential at -65 mV. 

The next generalization of these models would be to al- 
low both Gi and Ei, to vary independently, which would 
also change Er to varying degrees: for example, if G^ 



were very large, Er would be pulled towards E^,, while 
if Gi were very small, the value of i?i would have little 
effect on Er- Systematic investigations of the optimza- 
tions discussed here in such an extended model would 
be challenging due to the high computational costs of 
multidimensional parameter sweeps. 



C. Simulations 

All of our simulations were done using the NEU- 
RON/NMODL neuronal modeling language [l^ and 
auxiliary parameter-sweeping codes written in C and 
Python. We systematically varied G^ in Eq. ([3]) and oc- 
casionally other parameters, which are described in fur- 
ther detail in Section III, in order to determine how they 
influence the absolute and relative refractory periods of 
the action potential and the frequency of repetitive firing. 

Simulated action potentials were evoked in two dif- 
ferent ways. When studying the absolute and relative 
refractory periods, action potentials were evoked by 1 
A, 1 |as duration current injections into one end of the 
axon; we verified that these values were sufficiently brief 
and large so as to not influence the refractory periods. 
They are referred to throughout the text as "current 
spike-evoked" action potentials. The absolute refractory 
period. Tabs, was determined by finding the maximum 
time between successive current injections such that only 
one action potential resulted. We verified that for inter- 
injection times just above the absolute refractory period, 
two action potentials were produced and both propagated 
down the full length of the axon, as expected. The max- 
imum possible action potential firing frequency, i.e., the 
maximum frequency at which the axon can be driven 
with current spike inputs, is then the reciprocal of the 
absolute refractory period, fmax ■ 



fr, 



1 



bs 



(15) 



If the time interval between the input current spikes is 
greater than Tabs but less than a certain value Trei , called 
the relative refractory period, then while a second action 
potential is generated, it is generated in the wake of the 
first one, when the membrane and the voltage-gated ion 
channels have not yet returned to their resting states. 
The result is that the two action potentials interfere with 
each other: the velocity of the second is different from 
that of the first, and the distance between the action 
potentials - or, equivalently, the time between their peaks 
as measured at a single point along the axon - changes 
as they move down the axon. Information encoded in 
the distribution of intervals between action potentials can 
therefore be corrupted if they are too close together. If 
Ti is the time between the two current spikes and Tap 
is the time between the peaks of the two resulting action 
potentials as measured at some point further down the 
axon, then we can define the "interval shift" as AT = 
Tap - T,. For T, > Trei, Tap = T, and so AT = 0. 
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It is somewhat difBcult to determine Trei numerically. 
Action potentials in the Hodgkin-Huxely model have a 
phase of small, damped oscillations around the resting 
potential after the peak, which means that a second 
closely-following action potential can be either sped up 
or slowed down depending on which part of an oscillation 
it falls into. As a result, the interval shift AT itself oscil- 
lates around as a function of Ti (Fig. Inland [l^), and 
this causes T^ei as a function of G^, to have a jagged, dis- 
continuous appearance (Fig. IS]). Moreover, jitter noise 
in biological axons (caused by phenomena such as ion 
channel flicker) puts a nonzero lower bound on the tim- 
ing resolution of consecutive action potentials. 

Thus, it is both computationally easier and probably 
more biologically relevant to define a relative refractory 
period as a function of the maximum interval shift: for 
all values of Ti greater than Trei {^Tmax), by definition, 
I ATI = \Tap-T,\ < AT„^ax, where AT^ax > 0. The 
reciprocal of Trei (ATmax) gives the maximum frequency, 
frel {ATrnax), Sit which the axon can be driven with an 
interval shift no larger than AT^aa:: 



frel {ATmax) 



1 



Trei (AT„ 



(16) 



We determined Trei (ATmax) for values of ATmax ranging 
from 1 |xs (the numerical resolution of our simulations) 
up to 1 ms. 

To determine the repetitive firing frequency, /r, we 
simulated a constant (time-independent) current input, 
I DC, to one end of the axon. We then tested whether 
regular, repetitive firing resulted and, if so, measured 
the time between equivalent points on successive action 
potentials 20 ms after the beginning of the input current, 
by which time any initial transients had long since dis- 
appeared. The reciprocal of this time was then fr, which 
in general depended on the value of Idc as well as that 
of Gt. We refer to fr as the "repetitive firing frequency." 

We generally measured the time intervals between suc- 
cessive action potentials at a point 8 cm down the axon 
from the input stimuli. We verified that our optimiza- 
tion results were insensitive to the actual location of this 
point as long as it was outside a small region near the 
current injection site. 



III. RESULTS 



A. Resting Potential and Individual Action 
Potentials 



In addition to the effects of G^ on firing frequencies, 
which is the main focus of this study, also affects the 
shapes of individual action potentials and, in the case of 
CP channels, the overall resting potential of the axon. 
The resting potential Er is the voltage at which the sum 
of all the steady-state ionic currents is (Fig. [1}. It 
is therefore determined by the leak current as well as 
the small currents through the voltage-gated sodium and 
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FIG. 1: Despite its small magnitude, Gl has a significant 
effect on the resting potential of the axon, Er, defined as the 
potential at which there is no net transmembrane current. 
Here we plot Er as a function of Gl while keeping all other 
parameters fixed (this therefore corresponds to the case of 
Cl~ leak channels.) 



potassium channels, which are almost (but not entirely) 
closed at Er- Physically, the resting potential is deter- 
mined by the channel densities, the ionic concentration 
gradients, and the steady-state conformational configu- 
rations of the voltage-gated channels. 

Using the standard (and reasonably well-established) 
values of these other parameters, as Gl increases from 
0.05 to 3 mS/cm^, the resting potential increases (in the 
sense of getting less negative) by roughly 10 mV, with the 
sharpest rate of increase in the range below 1 mS/cm^. 
For Na'''/K~'' channels, there is no such dependence of Er 
on Gi because, as discussed above, we simultaneously 
varied the ratio of Na^ and leak conductances in 
order to keep Er at -65 mV. 

The effects of Gt on individual (current spike-evoked) 
action potentials differ between the two models. With 
Cl~ channels, higher values of G^ lead to smaller and 
narrower spikes that have a more pronounced post-peak 
oscillation. With Na~''/K~'' channels, by contrast, the 
heights and widths of the main action potential peaks are 
not significantly affected by Gl- In both models, smaller 
values of Gl cause the membrane potential to take longer 
to return to Er after an action potential (Fig. [5]); how- 
ever, the post-peak oscillations are larger for larger Gl 
in the Cl~ model, while they are smaller for larger Gl in 
the Na+/K+ model. 

We also investigated the effects of Gl on the metabolic 
energy consumption associated with action potentials, in 
view of other work [7, .8J suggesting that the overall scale 
of the three conductances in the squid giant axon is op- 
timized for the energy associated with action potential 
velocity. We did not find any such optimizations for Gl 
alone, but only a monotonic increase in metabolic energy 
consumption with increasing Gl- 
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FIG. 2: At lower values of the leak conductance, the axon 
takes longer to return to rest after a current spike-evoked 
action potential. Note that the resting potential is always 
made to be -65 mV for the Na^/K"*" case. 



B. Maximum Firing Frequency 

Figures [3] and |4] show some of our central results, the 
maximum firing frequency, fmax , calculated as a function 
of Gi. In the CP model, for values of much above 
the experimental range, fmax decreases by about 80 Hz 
for each 1 mS/cm^ increase in G^, (although the rate of 
decrease is slightly super-linear). 

However, as is just barely visible in Fig. [3] and much 
more evident in Fig. |4l the relationship between fmax 
and Gi in the CP model is not monotonic. For very low 
values of G^, fmax instead increases with G^, attaining a 
maximum value of about 560 Hz at 18.5 °C near G^, = 0.2 
mS / cm^ . Within the numerical limits of our simulation, 
the /moK-optimal value of G^ for the Cl~ model is about 
0.2 ± 0.06 mS/cm^, well within the range of experimen- 
tally measured values. 

In the Na^/K"*" model, however, fmax decreases mono- 
tonically by about 60 Hz for every 1 mS / cm^ increase in 
Gi. This may be due to the different way G^, affects 
the action potentials in this case (Fig. [2]). The optimal 
value of Gi in this model for maximum firing frequency 
is therefore 0. 

Since the rate coefficients (Eqs. (ISHTOl)) have a strong 
dependence on temperature, we repeated our calculations 
at 12.5 and 25 °C, representative of the range of tem- 
peratures the Loligo squid genus studied by Hodgkin and 
Huxley would normally experience. At both of these tem- 
peratures and for both of the models, we found qualita- 
tively similar behavior as at 18.5 °C. In the Cl~ model, 
the maximum values of fmax occur in or near the ex- 
perimental range of G^ (although the value of the max- 
imum fmax itself increases substantially with tempera- 
ture, from about 340 Hz at 12.5 °C to 848 Hz at 25 
°C). The optimal value of decreases with increas- 



FIG. 3: For Cl~ leak channels, the maximum firing frequency 
has a maximum value of about 560 Hz when Gl is around 0.2 
mS/cm^, in the range of experimentally measured values. By 
contrast, there is no local maximum for Na'^/K"'' leak chan- 
nels: fmax increases monotonically with decreasing Gl- Here 
and in subsequent figures, the two solid vertical lines show 
the range of Gl measured by Hodgkin and Huxley, while the 
dashed vertical line shows the mean value of their measure- 
ments. The Gl resolution in all simulations is 0.005 mS/cm^. 
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FIG. 4: The region around the Cl~ leak model fmax maxi- 
mum in Figure [3] is shown in greater resolution at two dif- 
ferent temperatures. The 18.5 °C case has its maximum 
at Gl ~ 0.2 ± 0.06. At 12.5 °C, the maximum fmax is at 
Gl ^ 0.27 ± 0.06 and has a value of about 340 Hz. The 
jagged shape of the curve and flatness near the minimum are 
due to the numerical limits of the simulation. 



ing temperature from Gl ~ 0.27 mS/cm'^ at 12.5 °C 
(Fig. 11) to Gl « 0.11 mS/cm2 at 25 °C. In both 
cases, the area around the maximum is fairly flat (al- 
though slightly less so for the 25 °C case) with a width 
of about 0.06 mS/cm^. While the optimal Gl is therefore 
not temperature- independent, the relative flatness of the 
maxima means that the maximum firing frequency for 
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values of Gl near, for example, 0.2 mS/cm^ are either at 
the maximum or within 2-3 Hz of it. 

In the Na^/K^ case, the relationship between and 
fmax at the other temperatures is still monotonically de- 
creasing with no local maximum. Additionally, no local 
maxima were observed when the magnitudes of the ac- 
tive sodium and potassium conductances, rather than the 
leak conductance, were varied. 



C. Repetitive Firing Frequency 

It is known ^2d\ from both theory and experiment 
that a constant current input to one end of a non-space- 
clamped axon can produce repetitive firing at a constant 
frequency, albeit only over a fairly narrow range of cur- 
rent. We investigated how the value of Gl affects this 
firing frequency, fr, which is qualitatively different than 
the firing produced by current spikes discussed above. 

For the Na'''/K~'' channel model, as with the maximum 
firing frequency, we found no evidence of an optimization 
of Gl at non-zero values for fr. The picture is consid- 
erably more complex, however, for the Cl~ model. Over 
a fairly large range of input current and temperature, fr 
attains its maximum value in or near the experimental 
range of G^- In Fig. [5l we show fr versus G^, for both 
models at typical values of temperature and input cur- 
rent. The Cl~ fr maximum in this case is about 208 Hz 
at Gl = 0.265 mS/cm^. 

However, the value of Idc is a second independent pa- 
rameter (assuming we hold all others fixed), and thus, in 
contrast to the maximum firing frequency produced by 
discrete current spikes (which is independent of their ac- 
tual size) , we must analyze the repetitive firing frequency 
as a function of both the leak conductance and the con- 
stant input current. In Fig. [BJ we show the repetitive 
firing frequency in the two-dimensional Gl-Idc parame- 
ter space. The bold curve shows the limits of the region 
where repetitive firing can occur. For combinations of 
Gl and Idc outside it, repetitive firing either does not 
occur or lasts for only a few spikes. The dashed curve 
within shows the /^-optimal value of Gl as a function of 
Idc', e.g., at Idc = 2.5 [a.A, the maximum fr of 218 Hz 
occurs when Gl is about 0.255 mS/cm^. If Gl is larger or 
smaller than this, fr decreases in a way similar to what 
is shown in Fig. [5l 

One evident feature of Fig. [5] is that no repetitive fir- 
ing at all is possible for values of Gl above about 0.6 
mS/cm^, regardless of the value of Idc- This maximum 
upper limit of Gl depends on temperature, as we will dis- 
cuss, but is always below about 1 mS/cm^ for tempera- 
tures above about 10 °C. Thus, Gl must be much smaller 
than the maximum voltage-gated sodium and potassium 
conductances in order for repetitive firing to occur. 

Another feature is that over about the lower half of 
the range of Idc where repetitive firing is possible, the 
optimal Gl is within the experimental range (0.13 to 0.5 
mS/cm^), and is everywhere below 0.27 mS/cm^. How- 
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FIG. 5: The repetitive firing frequency, fr, i.e. the action 
potential frequency caused by a constant current input (Idc) 
of 2.3 [lA at one end of the axon, is shown as a function of 
Gl- As before, no local maximum is evident for the Na"'"/K"'" 
leak channel model, and for it fr increases with decreasing 
Gl- However, the Cl~ model once again attains a maximum 
well above vanishing Gl- For this value of Idc and over a 
fairly large range in general, the maximum is very near the 
experimental value of Gl. 



ever, for values of Idc above about 3.8 |J.A, the maximum 
/,. obtains when Gl is 0. This includes the overall maxi- 
mum fr, 253 Hz, located at Idc — 4.08 [O-A. This value of 
fr is several tens of Hz above the maximum fr values at 
lower ranges where the optimal Gl is in the experimental 
range. 

In Figure [3 we show the combined results for an ap- 
proximately 10 °C range of temperature characteristic of 
what ocean-dwelling squid encounter over the course of 
a few months [2l|. At warmer temperatures, the range 
of Gl-Idc parameter space over which repetitive firing 
is possible decreases sharply. Conversely, the /,.-optimal 
Gl value as a function of Idc is relatively independent 
of temperature below the highest temperatures or above 
the lowest Idc values. Above about 4 |j.A, the optimal 
Gl is 0, with the highest absolute fr value also generally 
in this range, while for most of the range below, Gl is 
within the experimental limits. 



D. Relative Refractory Period 

As discussed above, the relative refractory period is 
more easily and probably more relevantly calculated with 
respect to a non-zero maximum interval shift ATmax - In 
Fig. [51 we show the corresponding maximum firing fre- 
quency, frei{A.Tmax), for a vcry low A.T,nax, 2 [is (inset), 
and then for two higher values, 200 and 300 \xs. The 
higher-AT,„aa; curves have optimal Gl values substan- 
tially outside the experimental range in the case of Cl~ 
channels, and only barely inside it for Na"*"/K''' channels. 
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FIG. 6: We show the behavior of the repetitive firing fre- 
quency in two-dimensional (Gl-Idc) parameter space. The 
soUd black curve shows the boundary of the region outside 
which a constant input current to one end of the axon does 
not produce sustained repetitive firing. The dashed curve in- 
side the boundary shows the optimal value of Gl for fr as 
a function of Idc- The numbers along the curve show the 
actual value of fr at various points. For the upper third of 
the Idc range inside the boundary, the highest value of fr is 
attained when Gl vanishes; for most of the rest, however, the 
optimal Gl is in the experimental range. 



For both models, the values of these optima are insensi- 
tive to ISTmax provided it is larger than about 150 |as. 
The Na+/K+ optima also appear to be relatively insen- 
sitive to other parameters such as temperature and axial 
resistivity. 

The low-ATmaa: curves look substantially different, and 
the reason for this is illustrated in Figure HI The relative 
refractory period for a given ISTmax is calculated by find- 
ing the largest value of the initial interval between current 
stimuli, Ti, such that |Ar| = ATmax! we then define this 
Ti value as Treii^Tmax), which can be visualized as the 
last point where the horizontal line representing l^Tmax 
intersects the Ti vs. |Ar| curve. For values of AT^ax 
larger than about 150 \is, this point always falls on the 
first, monotonically decreasing part of the curve below 
Ti ^ 5 ms, resulting in the smooth appearance of the 
higti-ATmax curves in Fig. [51 For low values of ATmax, 
however, the point of intersection falls within the region 
of secondary peaks due to the membrane potential os- 
cillations after the first action potential. The sizes and 
locations of the secondary peaks depend on Gl- As a 
result, it is possible for a peak to be just above AT^ax 
for one value of Gl and just below it for a second, nearby 
value of Gl, causing a sharp transition in the value of 
Trei{ATjnax) sls a fuuctiou of Gl- This is the case in 
the inset to Fig. [9] for Gl ~ 0.3 mS/cm^, the point 
of intersection is at about Ti = 16 ms, while for Gl = 
0.25 mS/cm^, the size of the peak there is slightly below 
ATrnax = 2 |js and therefore the intersection point jumps 



FIG. 7: We show the same information as in Figure [6] but 
for a range of temperature characteristic of what squid in the 
ocean encounter over periods of a few months. In general, 
higher temperatures decrease the size of the region of param- 
eter space in which repetitive firing occurs. The optimal Gl 
value for a given Idc is relatively temperature-independent 
except at high temperatures or low Idc values. 



down to the previous oscillation, at about Tj = 14.2 ms. 
We therefore see a sharp spike at Gl — 0.25 mS/cm^ in 
the inset to Fig. |8]for Cl~. Similar effects are seen in any 
curve of frei{ATmax) VS. Gl for small values of AT^ax- 



IV. DISCUSSION 

The two scenarios considered here for the stimula- 
tion of action potentials in the squid giant axon, delta 
function-like current spike inputs and unchanging con- 
stant current inputs, represent idealized extremes of the 
actual biology. The input current to the squid giant axon 
originates from postsynaptic glutamate-activated sodium 
channels at the squid giant synapse The frequency 
and duration with which these currents are evoked are ul- 
timately determined by the squid's sensory environment, 
e.g., whether it perceives any predators to be nearby. 
Thus, the actual current input to the axon when it is 
active is neither instantaneous nor constant, but is likely 
to be a time-varying function determined by the rate of 
synaptic bombardment and the kinetics of the synapse 
and postsynaptic sodium channels. Because of the rel- 
ative dearth of experimental data on the operation of 
the squid giant axon system in vivo (23j . it is difficult to 
accurately model this current. However, it might quali- 
tatively be expected that optimization results which hold 
for the two theoretical extremes would hold for the true 
time-dependent input current to the axon as well. 

Our results indicate that if the input stimulus to the 
axon can be regarded as a series of discrete, sharp pulses, 
and if the leak current is assumed to be chloride, then the 
maximum firing frequency of the axon is itself maximized 
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FIG. 8: The maximum firing frequency for a maximum al- 
lowed interval shift /\Tmax (defined in Eq. {16])) is shown for 
the two leak channel models for ATmax ~ 200 and 300 \is. 
For Cl~ channels, the frequency-optimal Gl is in the range 
of 1.2-1.3 mS/cm'^, well above the experimental range. For 
Na^ and channels, the optimal Gl is lower but still at 
the upper end of the experimental range. Inset: we show the 
frequencies for ATmax ~ 2 \is. At this as well as at other low 
values of the maximum interval shift, the frequency-optimal 
Gl values are well outside the experimental range in both the 
or (solid) and Na+/K+ (dashed) models. 



for values of leak conductance within the experimental 
range of values. For continuous input currents {Idc ), the 
same is generally true for at least half of the range of Idc 
and almost all the range of temperature over which repet- 
itive firing can occur. However, when one considers the 
full range of Gl-Idc parameter space, the repetitive fir- 
ing frequency is maximized overall at high values of Idc 
and Gl — Q- This raises the question of why, if the input 
current is constant (or, at least, has a timescale much 
longer than the repetitive firing frequency), a vanishing 
leak conductance and high input current level would not 
be preferable. 

One possible reason is metabolic energy consumption: 
higher I^c values would be associated with larger cur- 
rents through the Na^ channels in the synapse as well as 
more frequent action potentials in the axon itself. This 
implies that the metabolic energy cost of driving the axon 
at a frequency /r-, which can be quantified by the overall 
flux of Na^ ions into the membrane (all of which has to be 
subsequently pumped back out by the ATPase Na^/K^ 
exchanger in order to restore the resting concentration 
gradient), would be substantially higher at higher Idc 
and fr values. Of course, presumably it would not be 
optimal for fr itself to be too low. Hence, there may 
be an optimization involving both firing frequency and 
metabolic energy which would favor the lower end of the 
Idc range and the observed values of Gl- 

A simpler and perhaps more compelling reason is the 
requirement that the squid giant axon be able to func- 



FIG. 9: The absolute value of the interval shift, |Ar| = 
\Tap ~Ti\, is shown as a function of Ti for two nearby values 
of Gl- The secondary peaks after Ti = 5 ms are due to the 
post-peak oscillations of the first action potential, which can 
either speed up or slow down the second action potential. In- 
set: magnified view of the 10 ms < Ti < 20 ms region. The 
horizontal line shows |Ar| = ATmax = 2 ii.s; the largest Ti 
value at which it intersects each curve gives Trei{ATmax)- 



tion over a range of ocean temperature that can span 10 
or more degrees Celsius over the course of a few months 
[2l| . and several degrees over a single day [l^l- An im- 
portant consideration is that both Gl, which depends on 
the leak channel density on the axon, and IdCi which 
depends on the amount of neurotransmitter released per 
presynaptic action potential and the density of receptors 
on the postsynaptic membrane, are unlikely to be quickly 
changeable in response to a temperature change. (Even if 
the input current is not constant in time, its maximum or 
typical amplitude, of which we are taking Idc as a rough 
estimate, would depend on these properties and not be 
quickly changeable.) That is, for our purposes, we will 
assume that both Gl and Idc are essentially fixed. 

With these assumptions. Fig. [7] makes it clear that 
values of Idc in the vicinity of 3 |J.A are preferable to 
values much lower or higher in order for repetitive firing 
to be possible over the widest possible temperature range. 
For example, if Idc were 4 |j.A and Gl 0, which gives the 
maximum repetitive firing frequency at 18.5 °C, a modest 
temperature increase of only 1.5 °C would put Idc and 
Gl outside the region of parameter space where repetitive 
firing can occur. The axon would be rendered inoperable. 

Assuming that Idc niust be around 3 |J.A in order for 
the squid giant axon to function effectively at warmer 
temperatures means Gl should be at a value which in 
general gives the highest possible fr value at that Idc- 
We can see from Fig. [7|that over most of the temperature 
range, the optimal Gl values at Idc ~ 3 |a.A are clustered 
around 0.2 mS/cm^. Therefore, we may expect that Gl 
values in the experimental range are most optimal for 
repetitive firing given the typical temperature variations 
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the squid encounters. In sum, for Cl~ leak channels, 
biological values of G^, appear optimal for maximizing 
the maximum or repetitive firing frequency of the axon 
as determined by the Hodgkin-Huxley model. 

The same is not true, however, of Na~''/K"'" leak chan- 
nels, which show no such optimization at non-zero Gl 
values for the maximum or repetitive firing frequencies. 
They only appear to be superior in this regard to Cl~ 
channels when it comes to the relative refractory period 
as defined for a maximum allowed interval shift above 
about 150 [is: while both models evince non-zero op- 
tima for Gi, the one for Cl~ is far above the biological 
range of values, and the one for Na~''/K''' is only just 
within it. However, the calculation of relative refractory 
periods is problematic for lower values of the maximum 
interval shift. Moreover, the relevance of the relative re- 
fractory period, which characterizes the maximum firing 
frequency without information loss, to a peripheral axon 
like the squid giant axon is not clear. Thus, we consider 
our results for the firing frequency associated with the rel- 
ative refractory period to be much less compelling than 
those for the absolute and repetitive firing frequencies. 

V. CONCLUSIONS 

If one assumes that it is best for the squid's brain to 
be able to send two or more signals to its escape jet sys- 
tem with as little delay between them as possible, and 
also that the Hodgkin-Huxley model is a sufficiently ac- 
curate model of the biological squid giant axon, then the 
experimentally measured range of values for the squid 
giant axon leak conductance make far more sense for a 
chloride-like leak current than for one composed of sepa- 



rate sodium and potassium leak currents. The leak con- 
ductance appears to be optimal for for the firing rate of 
the axon, whether it be driven by discrete input current 
pulses or a by constant input current, if the leak current 
is assumed to be chloride or some combination of ions 
whose overall reversal potential is approximately Ea- If 
the leak current is instead assumed to be composed of 
separate sodium and potassium currents, then no such 
optimization is evident, though there is weak evidence of 
a partial optimization for relative refractory period. It 
should be remembered, though, that these results are all 
within the context of the Hodgkin-Huxley model. 

Of considerably more interest than the evolutionary 
neurobiology of Loligo is whether such an optimization of 
the leak conductance for firing rates exists or has taken 
place in mammalian central neurons. Due to the current 
lack of precise data on channel densities and kinetics in 
these much smaller and more morphologically complex 
cells, it is hard to address such questions in a rigorous 
way with modeling studies. Ncivcrtheless, whenever such 
data become available, it may be fruitful to examine in 
detail the role of the leak conductance on firing frequen- 
cies and information rates, as it appears to be a powerful 
mechanism for influencing these properties despite its de- 
ceptively small magnitude. 
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